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Abstract. Cepheids are the primary distance indicators 
for extragalactic astronomy and therefore are of very high 
astrophysical interest. Unfortunately, they are rare stars, 
situated very far from Earth. Though they are super- 
giants, their typical angular diameter is only a few mil- 
liarcseconds, making them very challenging targets even 
for long-baseline interferometers. We report observations 
that were obtained in the K' band (2-2.3 fj,m), on the 
Cepheid £ Geminorum with the FLUOR beam combiner, 
installed at the IOTA interferometer. The mean uniform 
disk angular diameter was measured to be 1.64 +0.14 - 
0.16 mas. Pulsational variations arc not detected at a 
significant statistical level, but future observations with 
longer baselines should allow a much better estimation 
of their amplitude. The distance to £ Gem is evaluated 
using Baade-Wesselink diameter determinations, giving a 
distance of 502 ± 88 pc. 

Key words: Stars: distances - Stars: individual: £ Gem 
- Stars: oscillations - Stars: Cepheids - Techniques: inter- 
ferometric 



1. Introduction 

Cepheids are very important stars in the history of astron- 
omy in general and cosmology in particular. The relation 
discovered by H. Leavitt (Pickering 1912) between the pe- 
riod of their light variation and their intrinsic luminosity 
makes them essential tools for the estimation of large as- 
tronomical distances. As they are intrinsically very bright 
stars, they can be detected in distant galaxies. The Hub- 
ble Space Telescope has successfully observed cepheids in 
M100, and measured a distance of 16.1 ±1.8 Mpc (Fer- 
rarese et al. 1996). The closest classical Cepheid (5 Cep) 
is situated at more than 300 pc. 

Their resolution by ground based instruments is a chal- 
lenge, as the angular diameter of nearby Cepheids is typi- 
cally smaller than 2 milliarcseconds (mas). In order to de- 
termine precisely the zero point of the period-luminosity 
relation, it is necessary to measure the angular diameter 
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of a number of Cepheids with high precision. The distance 
can then be derived from models of their intrinsic diam- 
eter with much better precision than via direct parallax 
measurement. The importance of Cepheid resolut ion b y 
interferometry is stressed in Sasselov & Karovska ( 1994 ). 

Until recently, the only available technique to measure 
the angular diameter of Cepheids was lunar occultation: 
it gave good results on ( Gem (Ridgway et al. 1982), with 
an angular size of 1.81 ± 0.31 mas in K. The progress of 
ground-based optical and infrared stellar interferometry 
allowed to measure a few Cepheid angular diameters in 
the last years (Mourard et al. 1997; Germain et al. 1999). 
We report in this paper our observations of £ Gem in the 
K' band with the IOTA interferometer equipped with the 
FLUOR beam combiner. 

In Sect. |^, we present the observations carried out at 
the IOTA interferometer that resulted in the C Gem data. 
In Sect. ^, we interpret these observations first in terms 
of mean angular diameter, and then we tentatively fit a 
variable diameter model to our data. In Sect. [|, using the 
extracted parameters, we can estimate by two different 
methods the distance to C Gem. We first use a Baade- 
Wesselink (BW) estimation of the intrinsic diameter of 
£ Gem. Another promising technique to estimate directly 
the distance to Cepheids is to use the angular amplitude 
of their pulsation coupled to radial velocimetry. Finally, 
Sect. H is dedicated to considerations for future observa- 
tions with the VLT Interferometer. 

2. Observations with FLUOR / IOTA 

2.1. Presentation of FLUOR and IOTA 

IOTA (Infrared and Optical Telescope Array) is located 
on top of Mount Hopkins, Arizona. It is a two telescope 
interferometer with 45 cm collecting apertures, and se- 
lectable baselines from 5 to 38 meters. It is operated both 
at visible and infrared wavelengths (Carleton et al. 1994; 
Traub 1998). A third telescope is currently being installed, 
together with a second set of delay lines. 

We have used the FLUOR (Fiber Linked Unit for Opti- 
cal Recombination) beam combiner. It is based on thin flu- 
oride glass single- mode waveguides (core diameter 6.5 jim) 
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for observations between 2 and 2.3 jim (K 7 band). A de- 
tailed description of FLUOR can be found in Coude du 
Foresto et al. ( 1998) ). The beam combiner accepts the light 
from two telescopes and produces four output signals: two 
photometric calibration signals (one for each beam) and 
two complementary interferometric outputs. Both spatial 
filtering of the input signals and beam combination are 
achieved simultaneously through the optical waveguides. 
Atmospheric corrugation of the incoming wavefront is fil- 
tered out at the injection in the fibers and converted into 
intensity fluctuations. The two photometric signals allow 
us to calibrate these variations continuously. 

This yields very accurate estimates of the squared 
modulus of the coherence factor /i 2 , which is linked to 
the object visibility V by the relationship 



Table 1. £ Gem and calibrator parameters 



/i 2 T 2 



(1) 



where T is the response of the system to a point source 
(interferometric efficiency). T is determined by bracket- 
ing the science target with observations of calibrator stars 
whose V is supposed to be known a priori. 

The current accuracy on visibility estimates with 
FLUOR is approximately 1 % for most sources and is as 
good as 0.3 % on bright objects. 

2.2. Q Gem and its calibrators 

£ Gem is the third brightest classical Cepheid in the north- 
ern sky [5 > -20 degrees) behind 8 Cep and rj Aql. More- 
over, £ Gem is the largest angular diameter Cepheid in the 
northern sky. This star is therefore particularly interesting 
for direct angular diameter measurements. 

During observations, the interferometric efficiency 
varies, which means that the science target observations 
have to be calibrated periodically using observations of a 
known, stable star. The calibration is then done by divid- 
ing the coherence factors obtained on the target by the 
interferometric efficiency values, derived from the calibra- 
tor observations. The choice of the calibrator is critical 
in the sense that any intrinsic variation in the calibrator 
visibility will contaminate the science target final data. 

The selected calibrators, HD 49968 and HD 62721, arc 
two bright K giants, located respectively 4.1 and 9.5 de- 
grees away from £ Gem (the parameters of all three stars 
are presented in Table Q). This proximity allows observa- 
tions of both stars in the same conditions, which limits the 
possible biases in the visibility calibration to a minimum. 
Practically, it also allows us to reduce the pointing time 
and to stay in the range of the delay line used to control 
the optical path difference. 

The spectral type of these stars make them very un- 
likely variable stars, and they are not reported variable in 
any catalogue of the SIMBAD database. Moreo ver, b oth 
calibrators have been selected by Cohen et al. ( 1999 ) as 
stable reference objects for infrared observations. These 



Name 


C Gem 


HD 49968 


HD 62721 




HR 2650 


HR 2533 


HR 3003 


rav 


3.62-4.18 


5.69 


4.88 


Mean mx 


1.98 


2.29 


1.32 


Spectral type 


F7Ib-G3Ib 


K5III 


K4III 


T off (K) 


5260-5780 


3800 


4000 


Log(Seff) 


1.5 


2.5 


2.5 


F B (M m ) 


13.11 


13.21 


13.21 


Parallax (mas) 


2.79 ± 0.81 


6.36 ± 0.92 


9.55 ± 0.83 


it 


0.282-0.241 


0.379 


0.370 


#ld/#ud 


1.021-1.018 


1.030 


1.029 


#ud (mas) 


~ 1.7 


1.87 ± 0.02 


2.93 ± 0.03 


Visibility 


~ 97 % 


~ 96 % 


~ 93 % 



authors have derived a limb-darkened (LD) angular diam- 
eter of 1.93 ± 0.02 mas for HD 49968 and 3.02 ± 0.03 
mas for HD 62721, based on high precision photometry 
at infrared wavelengths. Our model fitting is done with 
a uniform disk model, therefore it is necessary to convert 
these limb-darkened values to uniform disk (UD) equiva- 
lent diameters. 

The difference between the limb-darkened and uniform 
disk diameters is given by the relation (Hanbury-Brown et 
al. |1974D : 



7LD 



' l-u/3 
1 - 7u/15 



(2) 



with u = 0.379 and 0.370 the lin ear lim b darkening co- 
efficient in K (from Claret et al. 1995) respectively for 



HD 49968 and HD 62721. This coefficient depends on the 
effective temperature and gravity of the star (their val- 
ues are listed in Table |l|). This leads to a UD diameter of 
1.87 ± 0.02 mas for HD 49968, and of 2.93 ± 0.03 mas 
for HD 62721. This value differs from the limb-darkened 
diameter by 3 % , which is not negligible compared to 
the expected angular diameter variations of C Gem (about 
10 %). The C Gem limb darkening coefficient is also given 
for reference. It is variable as the temperature of the star 
changes during the pulsation. For simplicity, a mean value 
of 0ld/#ud =1.020 is retained for the LD values given in 
the rest of this paper. 

Other limb related effects, such as possible £ Gem limb 
brightening (proposed by Sasselov & Karovska 1994) are 
neglected. 

2.3. Visibility measurements 

Due to the very small expected angular diameter of C Gem, 
we used the longest available IOTA baseline, whose phys- 
ical length is 38 meters. The projected length on the sky 
during observations was about 37 meters. The precise spa- 
tial frequencies sampled during the observations are indi- 
cated in Table ||in cycles/arcsec. 



P. Kervella et al.: The angular diameter and distance of the Cepheid ( Gem 



3 



The observations were obtained in March 1999, De- 
cember 1999 and February 2000. £ Gem together with 
HD 49968 (during the 1999 campaigns) and HD 62721 
(for February 2000) were observed for sequences of twenty 
to thirty minutes alternatively. During each sequence, 100 
to 600 interferograms were acquired in a row. 

The data are acquired with FLUOR in scanning mode. 
A short stroke delay line (for the March 1999 cam- 
paign), later replac ed by a mirror mounted on a piezo 
stack (Ruilier |l999 installed in November 1999), sweeps 
through the zero optical path difference position, while a 
Nicmos-III based infrared camera measures the interfero- 
metric and photometric signals. For £ Gem, the fringe fre- 
quency was 100 Hz. Dark current sequences are recorded 
before and after each sequence and are used during the 
data reduction process for noise and signal calibration. 
Sequences are acquired every few seconds (March 1999) 
or twice per second (December 1999 and February 2000). 

We have used the method developed by Coude du 



Foresto et al. (1997) to reduce the data and to derive the 
coherence factor modulus of the interference fringes from 
each batch of observations. The interferometric efficiency 
of the instrument is computed from the measurements ob- 
tained on the calibrators. This process and the determi- 
nation of the object visibility are described in details in 
Perrin et al. (1998). The calibrated visibilities measured 
on ( Gem are presented in Table ^. A correction is applied 
to the calibrated visibilities to account for the difference in 
spectral type between the calibrators and ( Gem. The cor- 
rective factors, called shape factors (Fg), were computed 
for e ach star depending on its spectral type (Chagnon 
200C), and are listed in Table |l|. They are of the order 
of 13 /mi, and characterize the shape of the normalized 
spectral intensity distribution: 



Fr- 



B 2 (a)da 



(3) 



with B(o~) the spectral intensity distribution of the star 
normalized to unity: 



B(a)da = 1 



(4) 



The practical computation of Fb is achieved by integrat- 
ing the energy contained in the observation spectral band 
(K' in this pa per) o f standard stellar spect ra (tak en from 
Cohen et al. |l999| and Wallace & Hinkle |l997[ ). Before 
the integration, the spectrum is degraded to the actual 
spectral resolution of FLUOR (~ 50). At such a low reso- 
lution, the spectral lines disappear. This explains the slow 
variability with the spectral type (less than 1% between 
£ Gem and its calibrators). 

2-4- Data quality 

The data quality requirements for Cepheids observations 
are relatively higher than for other programs, as the angu- 



Table 2. Calibrated visibilities obtained with 
FLUOR/IOTA on ( Gem. 



Julian 


Phase 


Calibrated 


Sp.Freq. 


Scans 


Date - 




Visibility 


(cycles/ 




2451000 






arcsec) 




259.779 


0.3491 


0.9644 ± 0.0107 


83.76 


136 


262.722 


0.6391 


0.9777 ± 0.0108 


84.25 


157 


262.758 


0.6426 


0.9714 ± 0.0149 


83.78 


107 


595.838 


0.4581 


0.9932 ± 0.0085 


83.84 


316 


595.866 


0.4609 


0.9741 ± 0.0199 


83.78 


334 


602.734 


0.1375 


0.9700 ± 0.0116 


85.82 


396 


602.794 


0.1435 


0.9688 ± 0.0100 


84.20 


402 



lar sizes of these stars are very small, and their diameter 
variations are very subtle. Therefore, it is necessary to se- 
lect among the acquired data the best calibrated values to 
avoid the introduction of any bias in the star parameters 
evaluation. 

The data presented in this paper are the result of the 
selection of the highest quality measurements over the 
three observing runs. The selection was done according 
to the following rules: 

— £ Gem observations must be bracketed by a calibrator 
observation, to reduce the risk of a biased visibility 
evaluation. 

— the difference between the two successive transfer func- 
tion estimates (derived from the bracketing calibrator 
observations) shall not be larger than the sum of the 
half-error bars on those estimates: 

\T l+1 - T t \ < a{T t ) + a{T l+1 ) (5) 

— the x 2 between the calibrated visibility values from 
the two complementary interferometric channels must 
be lower than 3. A large discrepancy between the two 
channels would betray an instrumental or atmospheric 
perturbation. 

Only 1848 out of the 7459 interferograms acquired on 
C Gem (25 %) comply with all the quality conditions, and 
only those resulting visibility measurements are presented 
in Table [|. None of the data points obtained in December 
1999 was qualified, due to improper bracketing of (, Gem 
by calibrators. 

3. C Gem angular diameter 

3.1. Constant diameter model 

In this paragraph, we adopt a model of constant visibility 
for £ Gem. We then compute the mean angular diameter 
of this star over the observations. This is done by apply- 
ing a classical x 2 minimization algorithm with respect to 
the mean angular diameter only. The minimized quantity, 
relatively to the model diameter is 



x 2 = 



E 



(W0 



Vi„odol(0i)) 2 



(6) 
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Fig. 1. £ Gem radial velocity, taken from Bersier et al. 
( |19940 . The pulsation period is 10.150079 days. 

Table 3. Constant diameter model fit results 



Mean #ud 
Mean #ld 
Total x 2 
Reduced \ 2 



with 



1.66 +0.14 -0.16 mas 
1.69 +0.14 -0.16 mas 

6.08 

1.01 



V m 



jdcU 



7T e„ 



2Ji(«i) 



fi 



(7) 



(8) 



#mean is the mean uniform disk diameter in arcseconds 
and fi the spatial frequency in cycles/arcsec for measure- 
ment i, at the phase (f>i- The resulting mean uniform disk 
diameter is presented in Table 0. 

All the error bars in this paper are standard statisti- 
cal plus or minus one a (standard deviation) error bars, 
yielding a probability of 68 % for the extracted parame- 
ters to be in the error bars assuming no systematic bias is 
present. 



Nordgren et al. (2000), observing with the Navy Proto- 
type Optical Interferometer (NPOI), find a UD diameter 
value of 1.48 ± 0.08 mas (central wavelength Ao=0.735 
/zm). Independently, Lane et al. ( 2000| ) measured a mean 
UD diameter of 1.65 ± 0.3 mas with the Palomar Testbed 
Interferometer (PTI), in the H band (Ao=1.65 /im). Our 
mean UD diameter value is consistent with these two re- 
cent measurements. 

This estimate can also be compared with previous 
results from lunar occultation observations. Ridgway et 
al. ( 1982Q have measured the £ Gem diameter in the J 
and K bands. They have obtained uniform disk diame- 
ters of 1.81 ± 0.31 mas in K and 1.88 ± 0.86 mas in J. 
Chandrasekhar ( |1999| ) has measured a maximum value of 
1.9 ± 0.3 mas in the K band. Ashok et al. (1994) obtained 



a value of 1.6 ± 0.5 mas, also in K. These measurements 
and ours are compatible within the error bars though our 
value seems to be slightly smaller than the previous mea- 
surements. 
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Fig. 2. Diameter variations of £ Gem as integrated from 
the radial velocity data presented on Fig. [j], using a p- 
factor of 1.36. These variations are given relatively to the 
star radius at phase zero, when the maximum luminosity 
of the star occurs (for visible light). 



3.2. Variable diameter model 

3.2.1. Diameter variation curve from radial velocimetry 
data 

A more realistic model for £ Gem takes into account the 
angular diameter variations as they are measured by spec- 
troscopic radial velocimetry. In order to compute the curve 
of the radius variation of £ Gem, we integrated the high 
precision (er(v) < 0.5 km/s) radial velocity measurements 
obtained with the CORAVEL spectrograph by Bersier et 
al. (1994), presented in Fig. |l|. These measurements were 
phased using a period P = 10.150079 days and a reference 
epoch To = 2444932.736 taken from the same authors. The 



mean radial velocity of £ Gem (5.83 km/s) was subtracted 
from the radial velocity values before integration. 

The integration of the radial velocity curve requires 
to take into account the limb darkening of the star. The 
limb darkening gives a higher relative weight to the center 
of the disk than to the outer part. This means that the 
apparent radial velocity is lower than the real pulsational 
velocity value. The correction for this effect is included 
in a multiplicative term called the projection factor (or 
p- factor). 

The p-factor can be derived using models of Cepheid 
atmospheres (Sabbey et al. 1995 ) or using the measured 
widths and asymmetries of spectral lines. Its precise es- 
timation is still an open question. New interferometers 
under construction (see Sect, [s]) will enable direct mea- 
surements of the limb darkening on the closest Cepheids 
and hopefully will clarify this point. This coefficient is es- 
sential in the distance determination by the BW method 
(see Sect. 



11) 



The most widely accepted value for the 
p- factor of Cepheids is 1.36, and it has been used for our 
integration of the CORAVEL data. 

The resulting diameter curve is presented in Fig. ^| 
The shape of this curve is typical of an intermediate pe- 
riod Cepheid star, with a maximum at phase 0.3 and a 
minimum diameter at phase 0.7. For comparison, exam- 
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Table 4. Variable diameter model fit results 



Mean 8\jd 
Mean #ld 
Amplitude 
Total x 2 
Reduced \ 2 



1.64 +0.14 -0.16 mas 
1.67 +0.14 -0.16 mas 
0.1 ± 0.3 mas 

6.03 

1.21 



pies of typical radii curves for two Cepheids: U Sgr, P=6.7 
days and SZ Aql, P=17.1 days can be found in Laney et 
al. (1995). 

3.2.2. C Gem parameters 

We adjust a model with two parameters (the mean angu- 
lar diameter and the amplitude of the pulsation) to our 
visibility data. The shape of the expected angular diam- 
eter variation curve is shown in Fig. ^|. The fit is done 
in the visibility space. We use the same classical x 2 min- 
imization algorithm as for the constant diameter model, 
with two parameters. 

To extract the parameters of the £ Gem diameter vari- 
ations, the minimized expression is the same as in Eq. ^], 
but the model visibility is now: 



2Ji(z(&)) 



Knodel(</>i 

with 



Zfa) = * ( 0mean + ] /; 



AD U 



(9) 



(10) 



where: 



— #mean is the mean diameter of the star. This is the first 
parameter computed via the x 2 fit. 

— A is the pulsation amplitude, the second parameter 
adjusted by the x 2 fit. 

— xgk^ is the normalized diameter variation at phase 
4>i as given by the integration of the radial velocity 
curve. It is normalized to the total amplitude AZ? max . 

— fi is the spatial frequency in cycles/arcsec for measure- 
ment i. 

A plot of the result of the fit is presented in Fig. |3|. The 
parameters which minimize the x 2 are listed in Table |J. 



3.3. Discussion 

The reduced x 2 value for the variable model fit is 20 % 
larger than the value obtained for the constant diameter 
model. 

From this result we conclude that given the limited 
baseline of IOTA, the detection of the diameter variations 
is beyond the capabilities of our instrument. However, con- 
sidering the level of uncertainty on the individual points, 
both the constant and variable angular diameter models 
are compatible with our results. 
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Fig. 3. £ Gem variable angular diameter model fit. The 
curve corresponds to the parameters given in Table ^. 

4. Distance to £ Gem 

4-1- Using the Baade-Wesselink diameter estimate 

Following the suggestion from Sasselov & Karovska 
(1994), we compute the distance to £ Gem by combin- 
ing a mean radius from a complete B-W solution and our 
measurement of this star's angular diameter. 

The Baade-Wesselink (BW) method assumes that we 
can observe simultaneously the emitting surface of the 
star in flux, temperature variation and radial motion 



(Gautschy 1987 ). The flux and temperature variations give 
the ratio of the instantaneous radius of the star to the 
mean radius through the equation (assuming T e g is con- 
stant for simplicity): 



L(t 2 



10 



-0.2.(m(t 2 )-m(ti)) 



(11) 



The integration of the radial velocity curve gives the 
linear amplitude of the pulsation, making it possible to 
evaluate the radius itself. 

Numerous BW measurements have been carried out 
on £ Gem. A survey of eight BW diameter measurements 
obtained before 1982 can be found in Fernie (1984). They 
range from 64 to 73 Dq. More recently, Krockenberger 
et al. (1997), estimated the diameter of £ Gem to be 
69.1 +5.5 -4.8 Dq. Bersier et al. Ql997fl have derived a 
value of 89.5 ± 13.3 Dq. For coherence with the radial 
velocity measurements, this last value is retained in this 
paper's computations. 

Other methods, derived from BW, have also been used, 
such as 'CORS' by Ripepi et al. ( |1997| ), yielding a diameter 
of 73.5 Dq, or 'new CORS', giving 86.2 Dq. These authors 
have also experimented with variable values of the p-factor 
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on ( Gem during the pulsation, giving values of 80.1 Dq 
(CORS) or 93.7 D Q (new CORS). 

Knowing the mean uniform disk angular diame- 
ter of £ Gem and its intrinsic BW mean diameter, 



89.5 ± 13.3 Dq (from Bersier et al. 1997), it is now easy 
to compute its distance: 



9.305 D 



(12) 



and 



with d the distance in parsec, D the diameter in Dq, 
^ud the angular diameter in mas. 

The distance d to £ Gem can be derived knowing the 
constant diameter model 6>tjd = 1-66 +0.14 -0.16 mas. 
Eq. 12 gives d(( Gem) = 502 ± 88 pc, equivalent to a 



parallax of 1.99 ± 0.36 mas. The parallax obtained here 
for £ Gem is consistent with the parallax measured by 
Hipparcos : 2.79 ± 0.81 mas, though on the lower end of 
the error bars. 

It is interesting to note that the error in the derived 
distance is currently dominated by the uncertainty on the 
assumed BW mean diameter. Moreover, there is clearly a 
possible systematic error at least as large as the reported 
statistical error, considering the wide range of discrepant 
BW diameter values (64-94 Dq). The estimation of these 
systcmatics is beyond the scope of our paper and the error 
bars given here do not include them. 

4-2. Method for distance determination by the amplitude 
of the pulsation 

4.2.1. Principle 

It is also possible to estimate the distance to a pulsating 
star directly from the geometrical amplitude of its pul- 
sation. On one hand, we have the absolute amplitude of 
the pulsation by integrating the radial velocimetry curve. 
On the other hand, we can measure directly the angular 
diameter variation, though the amplitude determined in 
this paper is affected by a large uncertainty. By combin- 
ing these two values, one can derive the distance to the 
star via the equation: 



6.686.10" 9 AL> 



Ad 



I'D 



(13) 



with d in parsec, A#tjd m mas ano - in meters. 

Intuitively, we use what we know along the line of 
sight (radial velocimetry gives a linear amplitude value) 
and perpendicular to the line of sight (interferometry pro- 
vides an angular amplitude) to derive the distance. This 
assumes that £ G em is essentially pulsating radial ly (see 
Gautschy & Saio 1995 and Gautschy & Saio 1996 for an 
extensive review of star pulsations), but this hypothesis is 



nowadays widely accepted (Bono et al. 1999) 



4.2.2. Visible diameter curve vs. infrared interferometric 
observations 

While the radial velocity measurements are done mostly in 
the visible, the interferometric observations are often done 



at infrared wavelengths (Lane et al. |200fj, this paper). To 
assess the validity of the visible radius curve of the star to 
fit infrared wavelengths interferometric data, it is neces- 
sary to compare what is really observed by interferometry 
at 2.2 /Ltm to what is measured by radial velocimetry in 
the visible. 

The formation of the absorption lines used to measure 
the radial velocity through the Doppler shift happens at 
different depths in the stellar atmosphere depending on 
the wavelength: the longer the wavelength, the higher the 
line-forming region in the atmosphere. In the case of the 
infrared, the lines are formed above the visible ones. In 
the absence of shock waves, it is therefore expected to 
see a larger pulsation amplitude in the infrared than in 
the visible, as the infrared continuum-forming region is 
located above the visible photosphere. 



Regarding this question, Sasselov & Lester ( 1990 ) have 
obtained measurements at 1.08 and 1.6 /im of the radial 
velocity of several Cepheids, and find for £ Gem an appar- 
ent velocity difference of 2.7 km/s at phase 0.802, shortly 
after the minimum diameter. Confirming this trend, But- 
ler et al. ( |1997 ) have studied the cases of three bright 
Cepheids (S Cep, P=5 days; r\ Aql, P=7 days; X Cyg, 
P=16 days). They obtained radial velocity curves from 
spectral lines located in the near infrared (A = 0.8 fjxa) 
and in the visible (A = 0.5 /mi). For the two shorter pe- 
riod Cepheids, the near infrared and visible radial velocity 
curves differ by less than 1 km/s, but the longer period 
Cepheid X Cyg shows peak velocity differences of as much 
as 3 km/s. 

Assuming a radius for ( Gem of 89.5 Dq (1.25 10 11 m) 
taken from Bersier et al. (1997), a 3 km/s radial veloc- 
ity difference integrated over half the period would corre- 
spond to 2 % of the diameter of the star. The effect of this 
difference on the diameter curve is that the amplitude is 
increased by 20% (compared to the visible curve). Such a 
bias would also reduce accordingly the apparent distance 
estimation through the amplitude of the pulsation. 

Considering the fact that we do not detect the pul- 



sations significantly, this effect is neglected in Sect. 4.2.3 



However, when very high precision Cepheid diameter mea- 
surements will be possible with instru ments such as the 
VLT Interferometer (Glindemann et al. 2000 , see Sect. ||), 
it will be important to rely on radial velocimetry and inter- 
ferometric observations obtained at the same wavelengths 
for distance measurements. 

In the case of longer period, large amplitude Cepheids, 
the suspected presence of strong shock waves and velocity 



gradients (Butler et al. 1996) make it more difficult to esti- 



mate the true infrared pulsational velocity. Shorter period 
Cepheids are likely to be less affected, showing smoother, 
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more sinusoidal radial velocity curves than longer period 
stars. 



4.2.3. Application to £ Gem 

The amplitude of the pulsation as derived from the model 



fit presented in Sect. 3.2 is 0.1 ± 0.3 mas. The amplitude 
of the radius variation is 4.64 10 9 m (see Fig. ||), assuming 
that the atmosphere of the star is comoving (same behav- 
ior between visible and infrared wavelengths). We obtain 
d(C Gem) = 310 +oo -230 pc. 

This distance is compatible with the Hipparcos mea- 
surement (tt = 2.79 ± 0.81 mas, equivalent to a distance 
of 358 +147 -80 pc). Though this method is promising, 
the estimation of the distance to £ Gem by the pulsation 
amplitude is still beyond the present capabilities of our 
instrument. Moreover, due to the use of visible radial ve- 
locity curve s, our distance (as well as the one derived by 
Lane et al. |2000D may be bia sed by up to 20% towards 
smaller values (see Sect. 4.2.2) for discussion). 



5. Future observations with the VLT 
Interferometer 

In the next years, the VLT Interferometer (VLTI), 
CHARA and the Keck Interferometer will be the instru- 
ments providing the highest angular resolution. It is ex- 
pected that their hundreds of meters baselines will allow 
the resolution of the closest Cepheids. The VLT Interfer- 
ometer (up to 200 meter baseline, fringe tracking, large 
apertures with adaptive optics) will allow precise diam- 
eter measurements on a large number of Cepheids, cov- 
ering a wide range of periods and distances. Both ra- 
dial pulsations and non radial oscillations will be mea- 
surable through interferometric observations by using dif- 
ferent baseline orientations. 

Even assuming the same performances that were ob- 
tained with FLUOR (conservative 1 % visibility accuracy 
in 100 scans, without adaptive optics or fringe tracker), a 
study of the targets accessible to the VLTI gives a list of 
17 Cepheids for which the pulsation will be measurable at 
a 3 sigmas level per measurement. 



The VINCI instrument (Kervella et al. |2000D , function- 
ing on the same combination principle as FLUOR (single- 
mode infrared fibers) is expected to provide an even higher 
precision. 

One especially interesting target is 1 Car (HD84810). 
Its intrinsic diameter is evaluated between 180 Dq (Ripepi 
et al. 1997) and 195 D Q ( |1998[ ). Its relative proximity 
(Hipparcos parallax = 2.16 ± 0.81 mas) gives an expected 
angular diameter of more than 3.5 mas. On the 200 me- 
ters baseline of the VLTI, 1 Car will be fully resolved, and 
in the second lobe of its visibility function. It will there- 
fore be possible to look for features on the photosphere of 
this star. Since this star is a long period Cepheid (P=35.5 
days), it will be necessary to be particularly careful in 



the integration of the radial velocity curve to estimate the 
distance (see Sect. 4.2). 

C Gem will also be observable with the VLTI. Based on 
FLUOR observations, we can expect to measure the diam- 
eter of this star with a statistical precision of about 30 /ias 
per single measurement (100 interferograms, or about one 
minute of observation time) . Even considering calibration 
issues and the current poor knowledge of the p-factor, the 
final precision on the distance measurement should be bet- 
ter than 1 % for 300 observations spread over the pulsation 
period. 

It is also expected that direct limb darkening measure- 
ments will be possible on six Cepheids (for which visibility 
is less than 60 % on the 200 meter baseline). This will be 
achieved by observing the star at several baselines and 
computing the best fit with three parameters: limb dark- 
ening coefficient, mean angular diameter and amplitude of 
the pulsational variations. 

6. Conclusion 



We have reported the direct measurement of the angu- 
lar diameter of £ Geminorum at a relative precision of 
10 %, using long-baseline interferometry in the K' band. 
Our value is consistent with the previous determinations 
by interferometry (Nordgren et al. 2000, Lane et al. 2000) 
and lunar occultation (Ridgway et al. 1982). We have ten- 
tatively fit a two parameter model to our data, based on 
radial velocity measurements, in order to measure the am- 
plitude of the pulsation of this Cepheid, but the pulsation 
could not be detected on a significant statistical level. 

The distance to ( Gem was computed using BW lin- 
ear diameter estimate, and is consistent with the Hippar- 
cos parallax measurement. In case infrared interferometric 
measurements are used together with a visible light radial 
velocity curve, any direct distance estimation through the 
amplitude of the pulsation should take into account the 
difference between the pulsational motions of the visible 
and infrared photospheres of the star. Otherwise, an un- 
derestimation of the true distance could occur. 

The limb darkening factor is still not directly measured 
but estimated from a stellar atmosphere model. The use of 
large ground-based interferometers will allow to measure 
it directly. Once direct p-factor measurement and infrared 
radial velocimetry are available, the distance to the closest 
Cepheid stars will be determined in a fully observational 
way. 
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